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Identifying key players in collective dynamics remains a challenge in several research fields, from the 
efficient dissemination of ideas to drug target discovery in biomedical problems. The difficulty lies at several 
levels: how to single out the role of individual elements in such intermingled systems, or which is the best 
way to quantify their importance. Centrality measures describe a node's importance by its position in a 
network. The key issue obviated is that the contribution of a node to the collective behavior is not uniquely 
determined by the structure of the system but it is a result of the interplay between dynamics and network 
structure. We show that dynamical influence measures explicitly how strongly a node's dynamical state 
affects collective behavior. For critical spreading, dynamical influence targets nodes according to their 
spreading capabilities. For diffusive processes it quantifies how efficiently real systems may be controlled by 
manipulating a single node. 

Complex networks are a groundbreaking concept that is helping to understand the behavior of many 
chemical, biological, social and technical systems^'^. Network representations are particularly suitable 
for systems where heterogeneity dominates and are crucial for dynamics^, where a few nodes are usually 
considered as the most important. Oftentimes, node importance is correlated with centrality measures, local^'^ or 
globaP, which usually do not explicitly account for the dynamics as they are generally based on a purely 
topological perspective. However, dynamics is fundamental in assessing the impact of individual elements in 
global performance and in controllability problems^. Here, we show that dynamical influence is a centrality 
measure able to quantify how strongly a node's dynamical state affects the collective behavior of a system, taking 
explicitly into account the interplay between structure and dynamics in complex networks. We prove that it 
applies equally well to a variety of families of dynamical models, from spreading phenomena at the critical point to 
diffusive processes and and continuous-time dynamical system such as the Kuramoto model and the Roessler 
chaotic dynamics. 

Classical centrality measures in complex networks -like the degree or number of neighbors a node interacts 
with^'^, betweenness centrality^ counting the number of shortest paths through a certain node, eigenvector 
centrality^ based on the idea that relations with more influential neighbors confer greater importance, or the 
/c-shell decomposition^" that correlates with the outcome of supercritical spreading originating in specific nodes^^"^^ 
- rely only on topology, even if an underlying process can be indirectly associated in some cases. In contrast, the 
impact of individual elements in the global performance of the system inevitably depends on the specificities of the 
dynamics. Targeting individuals for vaccination strategies in epidemic processes is not the same as selecting 
electrical stimulation sites in the brain in order to suppress epileptic seizures. In this respect, a Laplacian-based 
centrality measure^^"^^, closely related to PageRank^^, has been proposed recently to assess the importance of 
complex network nodes in specific dynamical models. 

In this work, we provide a general and rigorous framework where dynamical influence is defined as a centrality 
measure both on directed and on undirected complex networks and applies to a variety of families of dynamical 
models, including epidemic spreading models like the susceptible-infected-removed (SIR), the susceptible- 
infected- susceptible (SIS), and the contact process, the Ising model, and diffusive processes like the voter model 
or phase coupled oscillators. In all cases, dynamical influence is calculated as the leading left eigenvector of a 
characteristic matrix that encodes the interplay between topology and dynamics. 



Results 

Defining dynamical influence. We focus on systems of AT time- dependent real variables x = (xi, 
coupled linear dynamics specified by a AT X AT real matrix M 



, x^) with 
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x = Mx. (1) 

A first classification of the dynamics is obtained by considering the 
largest eigenvalue fi^ax of For l^max < ^(t) converges to a null 
vector that represents a stable fixed point solution; for fi^^x > 
indefinite growth from almost all initial conditions is observed. 
Suppose that M is such that a non- degenerate /n^ax ^ 0 exists. 
Then, the scalar product (l)c = c ' x is a conserved quantity, where 
c is the left eigenvector of M for jUmax^ 

^=c-x{t) = [cM]-x{t) = 0. (2) 

The existence of the conserved quantity allows to calculate the final 
state in terms of the initial condition x(0) as 

/ N 1. / N C-X(O) . , 

x{(X)) : =\imt^oox{t) = e, (3) 

ce 

where e is a right eigenvector of M for jUmax- This equation implies 
that the projection of x{0) on c is all the system remembers at large 
times about the initial condition x(0). The coefficient q quantifies the 
extent to which the initial condition at node / affects the final state. 
Therefore, we call C/ the dynamical influence (DI) of element / on the 
dynamics under equation (1). 

One advantage of DI is that it is easily calculated without expensive 
numerical simulations. In fact, a simple way to calculate c furthers the 
understanding why this object quantifies the role of nodes in spread- 
ing dynamics. The power method (also called power iteration) 
approximates c by applying higher and higher powers of M to a 
uniform vector w^°^ = (1, 1, 1). For general exponent / E N, the 
/-th entry wf of 

w(^) = (l,...,l)M^ (4) 

is the number of all possible walks of length / departing from node / 
or, in other words, the number of ways an item can spread for / steps 
when originating at node /. At the first iteration this yields 

w^'^ = {l,...,l)M={d,A,...4n) (5) 

where di is the sum over the /-th row of M. When M is the adjacency 
matrix of a network, then wf^= di is the (out-) degree, the number of 

(2) 

(outgoing) connections of node /. For exponent 2, the /-th entry w] 
is the sum of the (out-) degrees of all neighbors of /. This is the same as 
the number of possibilities (walks) to depart from node / following 
two links. Now in the limit oo, the direction of the eigenvector c is 
approached by 

when the largest eigenvalue of M is non- degenerate and larger in 
magnitude than the other eigenvalues. Hence, the dynamic influence 
Ci of element / is its ability to serve as the origin of many arbitrarily 
long walks on the network. 

Epidemic spreading. Let us first apply these insights to critical 
phenomena like spreading processes In the SIR modeP°"^^, each 
node is either susceptible, infected or removed. An infected node / 
transfers the epidemic along each of its outgoing arcs independently 
with probability f^; node / itself relaxes to the removed state at unit 
rate. We study small perturbations to the stationary state with all 
nodes susceptible and approximate the dynamics by the linearization 

x=-x + ^A^x. (7) 

Here Xj{t) is the probability of node; to be infected at time t. The first 
term is the relaxation from the infected to the removed state at unit 
rate. The second term quantifies the transmission of the epidemic 
where the network enters by the transpose of its adjacency matrix A. 



Equation (7) can be rewritten as equation (1) with M = [^A^ - I, 
and J being the identity matrix. Matrix M has largest eigenvalue fi^ax 
= 0 when the spreading probability P is the inverse of the largest 
eigenvalue of A, that is P = fi^ = l/^max- We take again c as a left 
eigenvector of M at /Umax ^ 0 or, equivalently, a right eigenvector for 
maximum eigenvalue a^ax of A. Then the expected outbreak size 
from an initial infection described by the probability vector x{0) is 
proportional to c ' x{0). 

Now we ask how well c may forecast the actual SIR spread- 
ing dynamics, measured as the spreading efficiency (details in 
Methods) of node / that we define as the expected fraction of nodes 
reached by an epidemic outbreak initiated with node / infected {seed 
node), all others susceptible. Figure 1 shows that C/ is a good predictor 
of SIR spreading efficiency at critical parameter value P = in a 
small social network. Dynamical influence C/ outperforms the pre- 
dictions made by degree, shell index and betweenness centrality. 
Predictive power is quantified by the rank order correlation (see 
Methods). 

Figure 2 shows the predictive power of dynamical influence for 
spreading efficiency as a function of the infection probability in 
larger real-world networks and the Barabasi-Albert model. The 
results are as anticipated by the theory. Dynamical influence is a 
good predictor of spreading efficiency in the critical regime where 
P/Pc ~ 1- Predictions by dynamical influence outperform those by 
other quantities that are supposed to provide information about 
expected outbreak size in a broad interval of infection probabilities. 
This still holds for values of P that lead to average outbreak sizes of up 
to 10% of all nodes in the network, as indicated by the vertical dashed 
lines in Figure 2. 

The approximation w^^^ iox flnite length / in Eq. (4) is useful as a 
predictor of spreading efficiency as well. Even when the interaction 
network is not completely available, local information counting the 
number of walks of length I = 2 or I = 3 emanating from a node 
is enough to estimate dynamical influence. Figure SI in the 
Supplementary Information shows that the count of these short 
walks yields a prediction of spreading efficiency in the critical regime 
that is as good as dynamical influence itself The predictive power of 
these walk counts, too, reaches a maximum in the critical regime. 

For infection probabilities p far above or below the critical value P^, 
however, the degree di of a node / is a better predictor of spreading 
efficiency. In the subcritical regime, spreading is sparse and typically 
confined to the neighborhood of the seed node /, while in the super- 
critical regime, the epidemics rarely fails to spread to the whole 
system. In the critical regime in-between these extremes, infectious 
seeds are perturbations that trigger relaxation dynamics at all scales. 
This is reflected in a dynamics dominated by a marginal linear mode 
and a variety of possible final states. Dynamical predictions at crit- 
icality require then a global view of the network structure (and the 
final state is determined by the conservation law associated with the 
leading eigenvector c-removed). The scale-free distribution of epi- 
demic outbreaks in real populations^^'^^ is a sign of criticality and 
suggests that this regime is most relevant in practice. 

In order to check the robustness of the results we consider two 
modifications of the epidemic spreading dynamics. First, we study 
the SIR model with a stochastic rather than deterministic transition 
from the infected to the removed state. Specifically, the transition 
occurs with a probability fi independently for each infected node at 
each time step. Thus the time spent in the infected state (recovery 
time) has a geometric distribution with mean fi~\ This modification 
does not qualtitatively change the results of Figure 2 up to rescaling of 
P with fi. In fact, the curves of predictive power for different values of 
fi collapse when plotted as a function of the average outbreak size, see 
Figure S2. Second, prediction by dynamical influence may also be 
applied to the SIS model (see Materials and methods) yielding very 
similar results (Supplementary Information, Figure S3). The contact 
process^^ can also be considered, with A replaced by the stochastic 
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Figure 1 | SIR spreading efficiency compared to centrality measures in a social network. The network of Zachary's karate club^° has 77 edges connecting 
34 nodes, here ordered according to decreasing spreading efficiency. A monotonic decay of a centrality measure in the diagram indicates large predictive 
power for spreading efficiency. The rank order correlation of spreading efficiency is of 0.97 with dynamical influence, 0.86 with degree, 0.82 with shell 
index, and 0.79 with betweenness centrality. Indexing of nodes is the same as in ref.^°. In the network drawing, circles and squares represent the primary 
partitioning of the node set found by Girvan and Newman^\ Spreading efficiency has been estimated at j5 = = 0.15 performing 10^ independent runs of 
the SIR model per seed node. The largest eigenvalue of the network adjacency matrix is 6.65. 



adjacency matrix, the adjacency matrix after normalization such that 
each row sums up to 1. 

To facilitate intuitive understanding of the predictive role of cent- 
rality measures in spreading dynamics, let us consider small net- 
works. In each of the three cases in Fig. 3, a different subset of the 
measures yields the correct ranking by spreading efficiency at the 
critical point. The most efficient spreader is not necessarily the node 
with the largest degree. Being adjacent to several nodes with large 



degree may lead to large spreading efficiency despite a smaller degree, 
cf panel (a). This second order effect is reflected by dynamical influ- 
ence. When all degrees are equal as in panel (c), also dynamical 
influence and shell index are homogeneous. In this case, between- 
ness centrality captures the subtle effect of nodes having different 
positions in the network. We speculate that centrality measures 
based on unconstrained walks and shortest paths can do best in 
predicting spreading efficiency at the critical point. Then, a suitable 
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Figure 2 | Predictive power of different centrality measures for SIR spreading efficiency. Symbols are values of the rank order correlation coefficient of 
spreading efficiency with influence (squares), degree (triangles) and shell index (circles). The choice of the spreading parameter ^ controls the average 
outbreak size (horizontal axis), being the average number of nodes infected when choosing the seed node uniformly. The vertical dashed line indicates 
average outbreak size at the critical value of the spreading parameter ^ = j^^. The predictive power of betweenness centrality is below that of degree in all 
cases. The following networks have been used. E-mail interchanges between employees of a university^^, ~ 0.0482; unweighted neural circuitry of the 
roundworm C. elegans^^'^^ = 0.0654; snapshot of the Internet at Autonomous Systems level of Nov 08, 1997, see http://moat.nlanr.net, = 0.0315; a 
realization of the Barabasi-Albert (BA) model of scale-free networks^ with 1000 nodes and m = 2 edges added per node, — 0.0945. Other realizations of 
the BA model yield qualitatively the same result. For the BA model, shell index is not a predictor because its value ki = mis the same for all nodes. For the 
neural network, being directed, out-degree instead of degree is used as a predictor and for calculating the shell index. 
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Figure 3 | Comparison between spreading efficiency (diamonds), dynamical influence (squares) and betweenness centrality (stars) in small networks. 

In panel (a), the ranking of nodes with respect to spreading efficiency is rendered both by dynamical infiuence and betweenness. Note that the most 
efficient spreader is not a node with maximum degree but the node on the right connected to those maximum degree nodes. In the case of panel (b), the 
strongest spreaders are the nodes of maximum degree 3. However, the degree does not uniquely reveal the second strongest spreaders. Dynamical 
influence renders the full ranking of spreading efficiency. In panel (c), nodes are indistinguishable both by degree and dynamical infiuence. The small 
differences in spreading efficiency — note the scale on the axis — on this regular graph are rendered correctly by the betweenness centrality. The shell index 
is not usable as a node discriminator here. It takes value 1 on each node in panels (a) and (b) and the value 3 in panel (c). Spreading efficiency is calculated 
at the critical value P = j^^for each network, being 0.408 for (a), 0.463 for (b),and 0.333 for (c). For easier comparison, values of dynamical influence and 
betweenness centrality have been rescaled and shifted such that their mean and standard deviation are identical to that of the spreading efficiency in each 
network. 
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combination of dynamical influence with betweenness depending on 
network topology might be close to the optimal predictor. 

Nodes with large shell index are contained in highly connected 
neighborhoods that facilitate spreading. In many cases, the shell 
index may serve as a satisfactory predictor of spreading effi- 
ciency^^ Here, however, we find situations where its use for pre- 
diction is limited due to the degeneracy of the values shell index 
assumes. It has a constant value m across nodes of each network 
that can be built up by iterative attachment of a node with exactly 
m edges. This includes all trees {m = \) and networks from 
growth models such as the one by Barabasi and Albert^. A signifi- 
cant lack of resolution is also observed in real- world networks. 
Shell index assumes only few (—10) discrete values, cf. Figure S4 
in Supplementary Information. On the Internet graph, the same 
maximum shell index value is observed for the strongest spreader 
as well as nodes with spreading efficiency a factor of five below. 
Thus, even though the overall correlation between spreading effi- 
ciency and shell index is positive, lack of resolution limits the 
predictive power. Such limitations have been identified also in 
an empirical study of epidemic spreading in a social group^^, in 
a detailed comparison between SIR and SIS models^^ and in 
dynamics of rumour propagation^^ 

We remark that the most efficient spreaders are not necessarily the 
same as those targeted by efficient vaccination strategies in order to 
contain epidemics. At the network level, the aim of vaccination is to 
increase the epidemic threshold in order to render the spreading 
dynamics subcritical. The set of nodes by whose removal this shift 
of threshold is achieved^^ is different in general from the set of 
nodes with the largest dynamical influence. The Supplementary 
Information provides further results (Figure S5) and a brief discus- 
sion of vaccination. 



Ising model. The Ising model is a paradigmatic binary state model of 
critical phenomena. The Ising modeP°'^^ on a network^^ describes the 
dynamics of AT coupled spins S/ e { — 1; + 1} placed on the nodes. The 
zero temperature (T=0) version of the Ising model implements 
a majority rule for state updating. This is the same dynamics 
considered in threshold models of collective behavior for a 50% 
value of the threshold^^'^^, and its dynamics is also related to 
Schelling's model of urban segregation^^'^^; the finite temperature 
version has been considered in the context of strategic inter- 
actions^^. Finite temperature effects (noise), as considered here, 
are essential to escape from frozen configurations and to establish 
the robustness of transitions found in Ising-like models^^. Also 
in the theory of neural computation, Ising-like systems play an 
essential role^^. 

In the context of the Ising model, we define spreading efficiency of 
node / as the correlation between two measurements: the state of 
node / at time t and the magnetization (see Methods) of the whole 
system at a later time, formally 

m = {si{t)m{t + xm- (8) 

The parameter t measures the time lag between the two measure- 
ments. Figure 4 shows to which extent the ranking of nodes by Ising 
spreading efficiency is correlated with the ranking by various cent- 
rality measures. At the transition between order and disorder, Ising 
spreading efficiency has larger correlation with dynamical influence 
than with the other centrality measures. 

Diffusive processes: the voter model. Coming back to the general 
framework equation (1), there is a class of dynamical processes in 
networks in which the property of M having a zero maximum 
eigenvalue appears naturally without the need of adjusting any 





0.01 0.1 I 0.1 1 



average magnetization <lnil> average magnetization <lnil> 

Figure 4 | Predictive power of different centrality measures for Ising spreading efficiency at time lag r = 10 as a function of average absolute 
magnetization ( | /n | ). Symbols are values of the rank order correlation coefficient of spreading efficiency with influence (squares) , degree (triangles) and 
shell index (circles) . The vertical dashed lines indicate the value of ( I m I ) at the critical parameter value = f^c- Details on networks and the values of j^c are 
given in the caption of Figure 2. 
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parameter. This is the case of diffusive processes defined by equation 
(1) with M = —L and the Laplacian matrix entries 

N 

L^=-K^ + 3ijJ2Kik- (9) 

The zero eigenvalue of L is non- degenerate under mild 
assumptions^". For these processes our general analysis of equation 
(1) becomes exact. A prominent example of diffusive dynamics is the 
voter model^^ in which node / is in a spin state S/ e { — 1, + 1}. For this 
model, Xi stands for the ensemble average of spin /, Xi = (s^), and 
gives the rate at which node / copies the state of node j. Different 
definitions of the voter model dynamics provide clear examples 
of how the concept of dynamical influence takes into account the 
interplay between topology and dynamics: For link update dynamics 
in an undirected network, an ordered pair of nodes is chosen in 
each step and node / copies the state of node j. The rate matrix K 
becomes proportional to the transpose of the adjacency matrix A. As 
a consequence q = 1/N, the average magnetization X^^^^C/X/ is 
conserved, and all nodes have the same dynamical influence 
independently of the topological features of the network. In the 
more standard node update voter dynamics, at each step one node 
/ (having degree di) is selected at random and copies the state of one 
of its neighbors;, also selected at random. In this case Kij oc Ap/di, so 
that Kij is no longer a symmetric matrix, the conserved quantity is a 
weighted magnetization^^ and the dynamical influence of node / is 
proportional to its degree di. 

For diffusive processes, the system is driven towards a homogen- 
eous flnal state with x^'= = Xj{'^) for all / and;. Although x^ 
takes continuous values, each realization of the voter dynamics in a 
finite system eventually reaches a homogeneous absorbing state with 
either all nodes in the state + 1 or all in the state — 1. The influence C/ 
of a node weights the initial state of node / in the exit probability P+, 
that is, the probability to reach the absorbing configuration + 1 : P+ = 
((/>c + l)/2 = (x* + l)/2. When all nodes are equivalent (e.g., link 
update) X* is just the average of the initial values of the nodes, but 
otherwise (e.g. node update) x^ is given by a weighted average of the 
initial condition. The value C/ has an alternative interpretation as a 
stationary density of a random walk^^ 

Efficient driving of complex systems. The meaning of dynamical 
influence also manifests itself in the practical task of driving a system 
efficiently. In the context of the voter model, this task might be 
phrased in terms of the zealot problem^^. One considers a special 
directed network in which a given node (the zealot) does not copy the 
state of any of its neighbors. The question is the efficiency of the 
zealot in driving all other nodes to the zealot state. To show the broad 
applicability of the dynamical influence concept, we address this 
question of driving efficiency considering the problem of phase- 
coupled oscillators described by the Kuramoto model^^. Assuming 
all oscillators have the same intrinsic frequency w (without losing 
generality, we choose co = 0), the phase variable Xi of oscillator / 
advances as 

N 

Xi = co-\- Kij sin (xj — Xi). (10) 

;=i 

with a matrix K of non- negative coupling strengths. Around the 
synchronized state, phase differences are small. By approximating 
each sin-term with its argument, a linear homogeneous system as in 
equation (1) is recovered. 

We study a scenario with initially all oscillators / in phase X/(0) = 0. 
An additional node a with constant phase Xa = nil is added to the 
system and linked through an additional edge to a chosen node /. We 
measure the time T/ the system takes to reach the new homogeneous 
state with = nil for all nodes /. The dynamical evolution of these 
systems is illustrated by studying the motif in the inset of Fig. 5a. The 



global phase i/^(t) converges faster to the external forcing when the 
driving is applied to the nodes with higher infiuence, and the con- 
vergence of the different nodes depends on their relative network 
position in relation to the driver. In Fig. 5(b), we show the results on a 
directed network of phase oscillators connected as the network of 
regions in the macaque cortex^^. Dynamical influence has extremely 
high predictive power. The rank order correlation of driving effi- 
ciency with dynamical infiuence is 0.97, while 0.66 with degree ratio, 
— 0.14 with shell index and —0.09 with betweenness. Similar results 
are obtained on randomly grown directed networks and for coupled 
chaotic oscillators, see Figures S6 and S7 in Supplementary 
Information. These findings clearly show that dynamical influence 
is an excellent proxy to identify better targets for controlling global 
behavior, even in strongly non-linear dynamical systems. 

Discussion 

Dynamical influence is a centrality measure applicable to a wide 
range of dynamical processes on complex networks that takes into 
account the interplay between topology and dynamics. While the 
motivation and rigorous analysis of dynamical influence employ 
the context of linear systems, its practical use for understanding 
and controlling networked dynamics extends to several inherently 
non -linear systems. 

We have demonstrated that dynamical influence is applicable to 
stochastic equilibrium (Ising model) and nonequilibrium systems 
(epidemic and voter models) as well as deterministic state-continu- 
ous systems such as the Kuramoto model and the chaotic Roessler 
attractor. For critical epidemic spreading and the Ising model, 
dynamical influence is a good predictor of spreading capabilities. 
In the context of chaotic Boolean dynamics^^, a similar spectral cent- 
rality is highly correlated with a node's impact on the attractor 
reached^^. For diffusion, dynamical influence quantifies the impact 
of the dynamical states of single nodes on the asymptotic homogen- 
eous state. Beyond that, it proves to be a high-quality proxy for 
driving efficiency, uncovering which are the best target nodes in real 
networks to be forced in order to drive the system towards specific 
states. In a broader context, the identification of these targets has 
fundamental implications and practical applications on strategies 
with an interest in controlling collective behavior, from social influ- 
ence to biomedical responses. 

Methods 

Epidemic models. We simulate the SIR model of epidemic spreading in the time- 
discrete version. Transitions between the three states (S,I,R) are as follows. If node / is 
in the S (susceptible) state and has v infected (I) neighbors at time t, then node i 
remains susceptible with probability (1 — jS)", otherwise / is infected at time t + 1. If 
node i is in the infected state at time / then / is in the R (removed) state at time t+ 1. In 
the SIS model, at difference with SIR, a node infected at time t is susceptible again at 
time t + I. The probability of being removed in the SIR model does not enter in the 
linearized equation (7) because it appears only in a second order term in the equation 
for X. Therefore equation (7) gives the same linear description for the SIR and SIS 
models. 

The system is in an absorbing configuration if none of the nodes is infected. For 
both models, outbreak size is the number of nodes having been infected at least once 
before reaching an absorbing configuration. The spreading efficiency of node / is the 
average outbreak size when initiating the dynamics with node / infected and all others 
susceptible. 

Ising model. The spin values e { — 1, + 1} are updated asynchronously as follows. 
At each time step a node / ^ {!,..., AT} is drawn uniformly. The field 

h{t)=j2K,jS,{t) (11) 

;=i 

is calculated. The state of node i is flipped with probability 

min{exp[-jSsi(0/2i(0]'l}- (12) 

Flipping the state of node / means s,(f + 1) = —Si{t). Otherwise the state of node i 
remains unchanged. All other nodes j / retain their state, Sj{t + 1) = Sj{t). 

The parameter (inverse temperature) controls the order in the system. For large 
(small temperature), spins tend to align and there is long-range order seen as large 
clusters of nodes sharing the same spin value. For small (high temperature). 
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Figure 5 | Dynamical influence and driving in a system of phase-coupled oscillators, a, Adaptation of the global phase of the four-node system in the inset 
when driving is applied to one of the oscillator nodes (an additional oscillator coupled to the red circle node with fixed phase n/2. Adaptation is quick when 
driving at nodes 1 (o) and 4 (A), having high influence, and slow when driving at one of the other two nodes, having low influence, b. Driving 
efficiency (filled diamonds) and dynamical influence (shaded squares) for a system of phase oscillators coupled as the network of regions in the 
macaque cortex"^^. Driving efficiency of node i is measured as the time required to resynchronize with an additional input signal applied to a 
given node i. We say that the system has resynchronized when the global phase reaches^(t)>(l— )xa with a tolerance =10"^, where the global 
phase il/{t) is computed as the argument of the global order parameter z{t) = r exp[/i/^(^)] = ^^^.^ exp[zs2(t^]. For comparison, t^e degree ratio /Cout/^n of 
each node is also shown (open triangles). Since this is not an uncorrected network, node influence deviates significanfly from degree ratio. Each plotted 
quantity is rescaled by a factor to obtain an average value of 1. The empirical network serves as a testbed for prediction of driving efficiency. We do not aim to 
mimic real dynamics of the cortex. 



long-range order is absent. The magnetization 

m{t) = N-'Y,Si{t). 



(13) 



is used to quantif)^ the order of the system. Disordered systems have (| m|> ~ 0, while a 
finite positive value is obtained in ordered systems. 

Rank order correlation. For a vector x ^ IR", the rank of component / is given by 



r,{x) = \ + \{jiti\xj>Xi}\ + -\ {j it i\xj = x,}\ 



(14) 



The rank order correlation coefficient p{x,y) between two such vectors x and 7 is the 
Pearson correlation coefficient between the rank vectors r{x) and r{y). Thus p{x, y) 



takes values in [ — 1, 1] with p{x,y) -- 
increasing (decreasing) relation. 



+ 1 ( — 1) if and only if x and y are in a strictly 



Degree and degree ratio. The degree di of node / is the number of nodes / is connected 
to. In directed networks, in- and out-degree c/J" and d^^^ are distinguished. For the 
matrix averaging over all adjacency matrices of networks with fixed node degrees, = 
di is a left eigenvector for the largest eigenvalue. Likewise, the degree ratios c^f^*/ c^J" 
form a left eigenvector of the Laplacian matrix averaging over all networks with given 
node degrees^*^ '^''. 
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